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Low missing mass, single- and double diffraction dissociation is calculated for the LHC energies 
from a dual-Regge model, dominated by a Pomeron Regge pole exchange. The model reproduces the 
j /' j I rich resonance structure in the low missing mass Mx region. The diffractively excited states lie on 

the nucleon trajectory N* , appended by the isolated Roper resonance. Detailed predictions for the 
squared momentum transfer and missing mass dependence of the differential and integrated single- 
and double diffraction dissociation in the kinematical range of present and future LHC measurements 
are given. 

^ C"| The present work is a continuation and extension (e.g. with double diffraction) of a previous work 

O j [l[ using the dual Regge approach. 
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I. INTRODUCTION 



Measurements of single (SD), double(DD) and central(CD) diffraction dissociation is among the priorities of the 
. LHC research program. 

t-H 1 In the past, intensive studies of high-energy diffraction dissociation were performed at the Fermilab, on fixed 
deuteron target, and at the ISR, see 2] for an overview and relevant references. Fig. [3] shows representative curves 
of low-mass SD as measured at the Fermilab. One can see the rich resonance structure there, typical for low missing 
masses, often ignored by extrapolating whole region by a simple 1/M 2 dependence. When extrapolating (in energy), 
one should however bear in mind that, in the ISR region, secondary Reggeon contributions are still important (their 
relative contribution depends on momenta transfer considered), amounting to nearly 50% in the forward direction. 
At the LHC, however, their contribution in the nearly forward direction in negligible, i.e. less than the relevant error 



1 bars in the measured total cross section [3j. 

In most of the papers on the subject SD is calculated from the triple Regge limit of an inclusive reaction, as shown 
in Fig. [U 

In that limit, the double diffraction cross section can be written as 0, H[ 

d 2 a Gf 3 P 2 P(t) ( S \ 2 ^)- 2 ( Q H0)- Q p(t) 



dtdM 2 167t 2 Sq \s J \ s 



This approach has two shortcomings. The first one is that it leaves outside the small-M 2 resonance region. The 
second one is connected with the fact that whatever the Pomeron, the (partial) SD cross section overshoots the total 
one, thus obviously conflicting with unitarity. Various ways of resolving this deficiency are known from the literature, 
including the vanishing (decoupling) of the triple Pomeron coupling, but none of them can be considered completely 
satisfactory. 



We instead follow the idea put forward in paper Q and developed further in Ref. [lj, according to which the 
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FIG. 1: From SD to the triple Regge limit. 



Reggeon (here, the Pomeron) is similar to the photon and that the Reggeon-nucleon interaction is similar to deep- 
inelastic photon-nucleon scattering (DIS), with the replacement — Q 2 = q 2 — > t and s = W 2 — > M 2 . There is an 
obvious difference between the two: while the C parity of the photon is negative, it is positive for the Pomeron. We 
believe that while the dynamics is essentially invariant under the change of C, the difference between the two being 
accounted for by the proper choice of the parameters. Furthermore, while Jaroszewicz and Landshoff [5J, in their 
Pomeron-nucleon DIS structure function (SF) (or Pp total cross section) use the Regge asymptotic limit, we include 
also the low missing mass, resonance behavior. As is known, gauge invariance requires the DIS SF to vanish as Q 2 
(here, t) — > 0. This property is built in the SF, see Refs. [H, @|, and it has important consequences for the behavior 
of the resulting cross sections at low t (see Sec. IVII Appendix), not shared by the models based on the triple Regge 
limit, see d [|. 

It is evident that Regge factorization is essential in both approaches (triple Regge and the present one) . It is feasible 
when Regge singularities are isolated poles. While the pre-LHC data require the inclusion of secondary Reggeons, at 
the LHC we are in the fortunate situation of a single Pomeron exchange (Pomeron dominance) in the t channel in 
single and double diffraction (not necessarily so in central diffraction, to be treated elsewhere). Secondary Regge pole 
exchanges will appear however, in our dual- Regge treatment of Pp scattering (see below), not to be confused with 
the the t channel of pp. This new situation makes diffraction at the LHC unique in the sense that for the first time 
Regge-factorization is directly applicable. We make full use of it. 

The paper is organized as follows: Sec. [TT] contains a short overview on the pre-LHC measurements of (mainly) 
single and double diffraction dissociation and their theoretical interpretation. Sec. IIIII presents simple factorization 
relations that connecting DD with elastic and SD cross sections. Next, in Sec. IIVI we introduce the dual-Regge 
model the (non-linear) A^* Regge trajectory, the Roper resonance and summarize the main formulae used in our 
calculations. Fits to the data (figures an tables), presented in Sec. [V] are preceded by a brief introduction to our 
fitting strategy. In the Appendix (Sec. IVI|) we present some results of the main text but with the kinematical factor 
causing a dramatic turn down of the cross sections as t approaches 0. Sec. I VIII summarizes our conclusions. 



II. PRE-LHC MEASUREMENTS: CERN ISR AND SPS, FERMILAB. GENERALITIES 

Diffraction dissociation was predicted by theory (for a historical review and relevant references see @,[1,0|) and was 
intensively studied prior to the LHC on fixed targets (deuteron jets) at the Fermilab by a very successful Soviet-US 
collaboration Q and subsequently at the ISR, SPS and RHIC colliders, see 0, Q for more details. Below we briefly 
summarize the results and conclusions of those studies, mainly of single SD. 

1. First measurements concerned low missing masses near the threshold region. Resonance peaks e.g. in the reaction 
pp — >• X , depending on the relevant value of the momentum transfer t, were observed beyond the threshold. The 
origin and properties of theses peaks, still subject of debates, were revealed. The main deficiency of the missing mass 
method is the lack of any information on the decay properties of the produced particles. The missing information was 
obtained in bubble chamber experiments and in spark chambers, used also in exclusive channels at the ISR Collider. 

Apart from low missing masses, SD was found to persist also to high missing masses. Diffraction (coherence) 
imposes, however, an upper limit on the highest missing masses M x , roughly as £ = M\j s < 0.05. 

Beyond the resonance region, the smooth M x dependence of the cross section is approximated by d ^ M -i ~ M n , 
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hk-2. The value of n is related to the intercept of the trajectories exchanged in the proton- Reggeon scattering, as 
will be shown in Sec. IIV Al 

2. As expected, SD and DD are peaked in the forward direction, and the slope of the exponential peak of SD 
was found to be around 8 ~ 12 GeV~ 2 , varying with s, t and M 2 . Near the threshold, the slope is much (about 
twice) larger than that in elastic scattering, however near M x ^1.6 GeV it is already half of that of elastic pp. The 
correlation between the slope parameter and the mass of the excited state is a common feature of SD. 

A diffraction minimum around t rs 1 GeV 2 , similar to that in elastic hadron scattering, is expected also in SD (and 
DD). There are indications 10] of such a structure in pp — > p(nir + ) at rather small |t|, around \t\ ~ 0.2 0.3 GeV 2 
at yfs = 53 GeV, however its origin, fate and affinity to the dip-bump structure in elastic scattering is still a matter 
of debate Q. 

Diffraction is limited both in the missing mass (coherence), £ < 0.05 and in t ("soft" collisions). There is a 
transition region in t from "soft" to "hard" collisions, with a possible dip-bump structure between the two. To be 
sure, in our analysis we leave outside these interesting but controversial points, concentrating on the "first" cone with 
clear exponential behavior. 

Before the advent of the LHC, single diffraction dissociation was intensively studied in many different experiments: 
low energy ISR and SPS CERN experiments, low energy Fermilab (fixed deuteron target) experiments and high energy 
UA4, UA5, E710 and CDF experiments. All they cover the range 14 < ^/s < 1800 GeV, \t\ < 2 GeV 2 , and missing 

masses range from the threshold up to £ < 0.15. Here £ = — . The diffraction region is limited up to £ < 0.15. At 
£ ~ 0.15 non-diffraction contribution become sizable to diffractive one, and differential cross section become growing 
with £. The main results of these measurements and of their theoretical interpretation can be summarized as follows, 
for details see, e.g. @, [H|: 

1. Energy dependence. At energies below 30 GeV the integrated SD cross section rises with s according to the 
standard prescription of the the Regge-pole theory, however it slows down beyond. This effect was expected due to the 
familiar problem related to the violation of unitarity, namely that at high energies, implying the triple Pomeron limit, 
the DD cross section overshoot the total cross section, <tsd > &t(s)- Various means were suggested [l2| to remedy 
this deficiency, including decoupling (vanishing) of the triple Pomeron vertex. K. Goulianos instead renormalizes the 
standard Pomeron flux to meet the data, see [Tl|. Such a " renormalization" produces a break near y/s slowing down 
the rise of ctsd{s) in accord with the CDF data from the Tevatron, as shown in Fig. [21 taken from Ref. (TTj . 
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An alternative approach to resolve this crises (violation of unitarity) is possible within the dipole Pomeron approach 
Ref. [Hj]. 

2. t— dependence. SD cross section and the slope B(s,t, M 2 ) were measured in the range 0.01 < \t\ < 2. The 
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diffraction cone in SD essentially is exponential in t: a dip, similar to that in elastic scattering, is likely to appear 
(somewhere near t ~ 1 GcV 2 ). In previous work [lj, the cone of SD differential cross section shows a turn-down 
towards small-|i| due to the kinematical factor. Experimentally this turn-down was not yet confirmed. We ignore this 
effect in the main part of this work but will show it in the Appendix. This tiny effect is located in the kinematical 
region where Coulomb interaction is sizable, so it may hide the fine structure of the cone. However, as noticed in 
Ref. [HI], in SD (and DD), Coulomb interaction, at small squared momenta transfers, is suppressed compared to that 
in elastic diffraction pp scattering, allowing for a better determination of the strongly interacting part of the amplitude 
(in pp, at small \t\ this is possible only indirectly, by means of the Bethe-Heitler interference formula). 

3. M 2 dependence. Probably, this is the most delicate issue in the present studies (and diffraction in general). 
At the ISR, Fermilab and Tevatron, SD was measured in a wide span of the missing mass, starting from the inelastic 
threshold M 2 h = (m p + m^) 2 up to £ < 0.05 (or even to £ < 0.15 and higher), Fig.^b). 
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(a) (b) 
FIG. 3: Compilation of low-mass SD data form Fermilab experiments, see: @]. 



As shown in Figs. there is a rich resonance structure in the small M 2 region. In most of the papers on the 
subject, this resonance structure is ignored and replaced by a smooth function ~ M ~ 2 . Moreover, this simple power- 
like behavior is extended to the largest available missing masses. In Sees. IIVI and [V] we question this point on the 
following reasons: a) The low-M 2 , resonances introduce strong irregularities in the behavior of the resulting cross 
sections. LHC measurements, are able to probe low-M 2 region, and will be sensitive to these structures, b) The 
large- M behavior of the amplitude (cross sections) is another delicate point. Essentially, it is determined by the 
proton-Pomeron (pP) total cross section, proportional to the pP structure function, discussed in details in Sec. IIVI 
By duality, the averaged contribution from resonances sums up to produce high missing mass Regge behavior (M 2 )~™, 
where n is related to the intercept of the exchanged Reggeon and may be close (but not necessarily equal) to the 
above-mentioned empirical value ~ 1. 

In our calculations of the large- M behavior (Sees. IIVI and |V| we shall refer also to the ISR-Tcvatron-CDF data, 
remembering that for low energies 500 GeV), non-leading Regge exchanges are also important Q- 
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III. SIMPLE (AND APPROXIMATE) FACTORIZATION RELATIONS 



With the advent of the LHC, diffraction, elastic and inelastic scattering entered a new area, where it can be seen 
uncontaminated by non-diffraction events. In terms of the Regge-pole theory this means, that the scattering amplitude 
is completely determined by a Pomeron exchange, and in a simple-pole approximation, Regge factorization holds and 
it is of practical use! Remind that the Pomeron is not necessarily a simple pole: perturbative QCD suggests that the 
Pomeron is made of an infinite number of poles (useless in practice), and the unitarity condition requires corrections 
to the simple pole, whose calculation is far from unique. Instead a simple Pomeron pole approximation [16| is efficient 
in describing a variety of diffraction phenomena. 




FIG. 4: Diagrams for elastic scattering and diffraction dissociation (single, double and central). 



The D-L elastic scattering amplitude is simply [16| 

a( S , t) = mm^s/sor^- 1 + a r ( S , *), a) 

where is the signature factor, and a(t) is the (linear) Pomeron trajectory. The signature factor can be written 
as £(i) = e - "/ 2 , however it is irrelevant here, since below we shall use only cross sections (squared modules of the 
amplitude), where it reduces to unity. The residue is chosen to be a simple exponential, f3(t) = e bpt . "Minus one" in 
the propagator term (s/so) Qp< -'- l_1 of (JT|) correspond for normalization <tt(s) — ImA(s, t = 0). The scale parameter Sq 
is not fixed by the Regge-pole theory: it can be fitted do the data or fixed to a "plausible" value of a hadronic mass, or 
to the inverse "string tension" (inverse of the Pomeron slope), so = 1/a'. The second term in Eq. ([T]), corresponding 
to sub-leading Reggeons, has the same functional form as the first one (that of the Pomeron), just the values of the 
parameters differ. We ignore this term for reason mentioned above. 

Fig. 2] shows the simplest configurations of Regge-pole diagrams for elastic, single- and double diffraction dissocia- 
tion, as well as central diffraction dissociation (CD). In this paper we consider only SD and DD. 

Factorization of the Regge residue j3{t) and the "propagator" (s/so) ap ^~ 1 is a basic property of the theory. As 
mentioned, at the LHC for the first time, we have the opportunity to test directly Regge-factorization in diffraction, 
since the scattering amplitude here is dominated by a simple Pomeron-pole exchange, identical in elastic and inelastic 
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diffraction. Simpie factorization relations between elastic (-ff 1 ), single (-ff 1 ) and double ( dtdM°dM' 2 ) are known 
from the literature Really, by writing the scattering amplitude as product of the vertices, elastic / and inelastic 
F, multiplied by the (universal) propagator (Pomeron exchange), f 2 s a fFs a , F 2 s a for elastic scattering, SD and 
DD, respectively, one gets 

d 3 a D D _ d 2 a S Di d 2 a S D2 , d(7 el 
dtdM 2 dM 2 ~ dtdM 2 dtdM 2 1 dt ■ U 



Assuming e Bt as a t— dependence for both SD and elastic scattering, integration over t yields: 

d 3 a DD _ d 2 a S Di d 2 a SD2 ; , , 

dMfdM 2 ~ dM 2 dM 2 /<7el - [ j 

where k = r 2 /{2r - 1), r = B SD /B el . 

For pp interactions at the ISR, r — 2/3 and hence k = 4/3. Taking the value r = 2/3, consistent with the 
experimental results at Fermilab and ISR, one obtains <tdd — r Notice that for r — 1/2, k — > oo. Thus k is very 
sensitive to the ratio r, which shows that direct measurements of the slopes at the LHC are important. Interestingly, 
relation ([3]) can be used in different ways, e.g. to cross check any among the four inputs. 

To summarize this discussion, we emphasize the important role of the ratio between the inelastic and elastic slope, 
which at the LHC is close to its critical value Bsd/B c i = 0.5 (it cannot go below!), which means a very sensitive 
correlation between these two quantities. The right balance may require a correlated study of the two by keeping the 
ratio above 0.5. This constrain may guide future experiments on elastic and inelastic diffraction. For detailed details 
see also [13, El- 



IV. MODEL FOR SINGLE AND DOUBLE DIFFRACTION DISSOCIATION 



The model relies on the following premises: 

1. Regge factorization is feasible since, as stressed repeatedly, at the LHC energies in the region of \t\ < 1 GeV 2 , 
which is typical for diffraction, the contribution from secondary Reggeons is negligible, and, for a single Pomeron 
term, factorization ^ is exact. Due to factorization, the relevant expressions for the cross sections (elastic, SD, DD) 
have simple forms ©, (ITU1) . (fTTj) . Such relations are known from the literature, see e.g. 0, [lj], [H| and references 
therein. 

2. The inelastic pPX vertex receives special care. 

Having justified and accepted the factorized form of the scattering amplitude, the main object of our study is now 
the inelastic proton-Pomeron vertex or transition amplitude. As argued in Refs. P, @, it can be treated as the 
proton structure function (SF), probed by the Pomeron, and proportional to the Pomeron- proton total cross section, 
al^ p (M 2 ,t), with the norm al^ p (M 2 ,t) = ImA(M 2 ,t), in analogy with the proton SF probed by a photon (in ep 
scattering e.g. at HERA or JLab). 

,W 2 (Mlt) = F 2 ( X ,t) = — - ^ 2 a/H))3/2 W(M»,t), 

where a is the fine structure constant, v = Mx 2 J^ p - , and x — 2 m u 1S * nc Rjorken variable. 

The only difference is that the Pomeron's (positive) C parity is opposite to that of the photon. This difference 
is evident in the values of the parameters but is unlikely to affect the functional form of the SF itself, for which we 
choose its high-M 2 (low Bjorken x) behavior. Notice that the the total energy in this subprocess, the analogy of 
s = W 2 in DIS, here is M 2 and t here replaces q 2 = —Q 2 of DIS. Notice that gauge invariance requires that the SF 
vanishes towards Q 2 — > (here, t), resulting in the dramatic vanishing of the SD and DD differential cross section 
towards t = 0. How fast does the SF (and relevant cross sections) recover from t = a priori is not known. 
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FIG. 5: Virtual photon + proton — > M% transition. 

Furthermore, according to the ideas of two-component duality, see e.g. 19], the cross sections of any process, 
including that of pP — > X, is a sum of a non-diffraction component, in which resonances sum up in high-energy 
(here: mass M 2 plays the role of energy s) Regge exchanges and the smooth background (below the resonances), dual 
to the Pomeron exchange. The dual properties of diffraction dissociation can be quantified also by finite mass sum 
rules, see [2|. In short: the high-mass behavior of the pP — > X cross section is a sum of a decreasing term going 
like ~ jj^,m ss 2 and a "Pomeron exchange" increasing slowly with mass. All this has little affect on the low-mass 
behavior at the LHC, however normalization implies calculation of cross sections integrated over all physical values 
of M 2 , i.e. until M 2 < 0.05s. 

3. The background is a delicate issue. In the reactions (SD, DD) under consideration there are two sources 
of the background. The first is that related to the t channel exchange in Fig [3](b) and it can be accounted for by 
rescaling the parameter sq in the denominator of the Pomeron propagator. In any case, at high energies, those of 
the LHC, this background is included automatically in the Pomeron. The second component of background comes 
from the subprocesses pP — > X. Its high-mass behavior is not known experimentally and it can be only conjecture on 
the bases of the known energy dependence of the typical meson-baryon processes appended by the ideas of duality. 
The conclusion is that the Pp total cross section at high energies (here: missing masses M) has two components: a 

decreasing one, dual to direct-channel resonances and going as af£ ~ ^2r( s ')°' r ^~ 1 = J2r (M 2 ) " 1 ^ 1 , where R 
are non-leading Reggeons, and a slowly rising Pomeron term producing ~ M 2 ' 08 (l6| . 



A. Duality 



Any meson-baryon total cross section (or scattering amplitude) is a sum of two contributions: diffractive and non- 
diffractive. By the concept of two-component duality 20], the diffractive component, the smooth background at low 
energies (here: missing masses) is dual to a Pomeron i-channel exchange at high energies, while the non-diffractive 
component contains direct channel resonances, dual to high-energy t— channel (sub-leading) Reggeon exchanges, as 
shown in Fig. [51 




Veneziano duality 



FIG. 6: Connection, through unitarity (generalized optical theorem) and Veneziano-duality, between the inelastic form factor 
and sum of direct-channel resonances. 



According to our present knowledge about two-body hadronic reactions, two distinct classes of reaction mechanisms 
exist. 
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The first one includes the formation of resonances in the s— channel and the exchange of particles, resonances, or 
Regge trajectories in the t— channel. The low-energy resonance behavior and the high-energy Regge asymptotics are 
related by duality, which at Born level, or, alternatively, for tree diagrams, mathematically can be formalized in the 
Veneziano model, which is a combination of Eulcr Beta- functions [21j . 

The second class of mechanisms does not exhibit resonances at low energies and its high-energy behavior is governed 
by the exchange of a vacuum Regge trajectory, the Pomeron, with an intercept equal to or slightly greater than 
one. Harari and Rosner [2(| hypothesized that the low-energy non-resonating background is dual to the high-energy 
Pomeron exchange, or diffraction. In other words, the low energy background should extrapolate to high-energy 
diffraction in the same way as the sum of narrow resonances sum up to produce Regge behaviour. However, contrary 
to the case of narrow resonances, the Veneziano amplitude, by construction, cannot be applied to (infinitely) broad 
resonances. This becomes possible in a generalization of narrow resonance dual models called dual amplitudes with 
Mandelstam analyticity (DAM A), allowing for (infinitely) broad resonances, see |22| and references therein, or the 
background. 

In the resonance region, roughly 1 < M < 4 GeV, the non-diffractive component of the amplitude is adequately 
described by a "reggeized Breit-Wigner" term Eq. following from the low-energy decomposition of a dual ampli- 
tude, with a direct-channel meson-baryon (Pomeron-proton, in our case) trajectory (see Fig. [7] in Sec. HVBp . with 
relevant nucleon resonances lying on it, appended by a Roper resonance (see Eq. ([5])). 



XmA{a + b -> c + d) = 


R 


Pomeron 


s— channel 




Non-resonant background 


t— channel 




Pomeron (I = S = B = 0; G = +1) 


High energy dependence 


s"" 1 , a < 1 


s a ~\ a > 1 



TABLE I: Two-component duality 



By duality, a proper sum of direct channel resonances produces smooth Regge behavior, and, to avoid "double 
counting" , one should not add the two. Actually, this is true only for an infinite number of resonance poles. For 
technical reasons, we include only a finite number of resonance poles, moreover, apart from the "regular" contribution 
of the nucleon resonances, lying on the N* trajectory, the Roper resonance is also included (see Sec. HVCp . The 
spectroscopic status of this resonance is disputable. It has no place on the known baryonic trajectories, but its hight 
exceeds that of the neighboring, next most important TV* (1680). In view of the "truncated" series of resonance poles, 
we do not expect that it will reproduce correctly the high-energy Regge behavior, therefore we add it to the total 
cross section in the form of an "effective" Regge pole contribution. 

The second, diffractive component, is essentially the contribution from a Pomeron pole exchange a^, p ~ (M 2 ) a ( ) — 1 , 
with a(0) ss 1.08 [lj|, and as shown by Donnachie and Landshoff [l6| this term also can give some contribution to 
the low-energy (here, missing mass) fiat background. 

To summarize this discussion, Regge-pole exchanges take place at two distinct parts of the diagrams shown in Fig. |3 
in the t channel, where, at the LHC, only the Pomeron contributes, and in the inelastic form factor, sub-diagram 
shown in Fig. [51 where, depending of the value of the missing mass, both the Pomeron and Reggeons are equally 
important. At low missing masses, the direct-channel proton trajectory N*, dominates, replaced by an effective 
Reggeon exchange at high masses, appended by a Pomeron. The Roper resonance (with its controversial status) in 
the direct channel stay apart. Although we concentrate on the low missing mass region, the behavior of the cross 
sections at high masses are important in the calculation of the cross sections integrated in M 2 . We remind that we 
impose a limit on diffraction events to be about £ < 0.05 or M < 200 GeV (as it used at some experiments). 



B. Resonances in the Pp system; the TV* trajectory 

The Pp total cross section at low missing masses is dominated by nucleon resonances. In the dual-Regge approach 
[l|, the relevant cross section is a "Breit-Wigner" sum Eq. (jj), in which the direct-channel trajectory is that of N* . 
An explicit model of the nonlinear N* trajectory was elaborated in paper [23[ and used in [l[. 
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cr£ p {Ml ,t)=Xm A(M 2 X ,t)=Xm( 

Vn=0,l,. 



af(t) 



2(n+l) 



2n + 0.5 -a N .{M%) 



(4) 



The Pomeron-proton channel, Pp —> Ml couples to the proton trajectory, with the I(J P ) resonances: 
1/2(5/2+), F15, m = 1680 MeV, T = 130 MeV; 1/2(9/2+), H w , m = 2200 MeV, T = 400 MeV; and 
1/2(13/2+), K hl3 , m = 2700 MeV, T = 350 MeV. The status of the first two is firmly established [H, while 
the third one, iV*(2700), is less certain, with its width varying between 350 ± 50 and 900 ± 150 MeV . Still, with the 
stable proton included, we have a fairly rich trajectory, a(M 2 ), whose real part is shown in Fig. [3 

Despite the seemingly linear form of the trajectory, it is not that: the trajectory must contain an imaginary part 
corresponding to the finite widths of the resonances on it. The non-trivial problem of combining the nearly linear and 
real function with its imaginary part was solved in Ref. [23j j by means of dispersion relations. 
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FIG. 7: The real part of the proton Regge trajectory. The 
dashed line corresponds to the result of a linear fit, the solid 
line is the fit from I23I. 



FIG. 8: The widths of the resonances, F : 



Xm q(M- i ) 
MTZea'(M 2 ) ' 

ap pea ring on the proton trajectory, calculated and fitted 
in EH. 



We use the explicit form of the trajectory derived in Ref. 



ensuring correct behaviour of both its real and 



imaginary parts. The imaginary part of the trajectory can be written in the following way: 



Xm a(s) = s s c n 



[s - s n ) , 



(5) 



where A n = IZe a(s n ). Eq. ([5]) has the correct threshold behaviour, while analyticity requires that S < 1. The 
boundedness of a(s) for s oo follows from the condition that the amplitude, in the Regge form, should have no 
essential singularity at infinity in the cut plane. 



The real part of the proton trajectory is given by 

TZea(s) = a(0) 



(6) 
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where 

i / \ r(l- 5)T(Xn +1) ( , „ s 

I (A„ - + 2Js„ V s„ 

cot[7r(l - 6)}- 



TTS 6 - 1 



r(-<5)r(A„ + i)s 



s r(A„ - s + 



^2F 1 (6 - A„, 1; 5+ 1; ^) | 0(s - s n ) 



The proton trajectory, also called N+ trajectory [23|, contains the baryons N(939) ± + , N(1680) | + , N(2220) 9 + 



and N(2700) " [24|. In the fit, the input data are the masses and widths of the resonances. The quantities to be 



determined are the parameters c n , S and the thresholds s n . Following |23j we set n = 1,2, a; and si = (m n + m^) 2 = 
1.16 GeV 2 , s 2 = 2.44 GeV 2 and s x = 11.7 GeV 2 . 

Other parameters of the trajectory, obtained in the fit, are summarized below: a(0) = —0.41, 6 = —0.46 ± 0.07, 
c\ = 0.51 ± 0.08, C2 = 4.0 ± 0.8 and c x — (4.6 ± 1.7) • 10 3 . Taking the central values of these parameters we obtain 
the following values for the A's: A x = 0.846, A 2 = 2.082, \ x = 11.177. 

The elastic contribution, is separated from SD and DD by a gap extending from the proton mass m p to the first 
threshold, at m p + m T , and it should be treated separately. 



Thus, we obtain: 



i mM MU) = .EWlf"- (2n+a5 _ affjgffi (Im „ ■ P) 



C. The Roper resonance 



Apart from the well established protonic trajectory with a sequence of four particles, there is a prominent single 
resonance / = 1/2, J = 1/2+ with mass 1440 MeV, known as the Roper resonance [24j . It is wide, the width being 
nearly one fourth of its mass, its spectroscopic status being disputable. There is no room for the Roper resonance on 
the proton trajectory of Sec. IIVB1 although it could still be a member of protons daughter trajectory. Waiting for a 
future better understanding of Roper's status, here we present the contribution to SD cross section of a single Roper 
resonance, calculated from a simple Breit-Wigner formula: 



lmA Roper (M x ,t) = b ■ — -WZTr~ Joy*' ( 8 ) 

\Mx ~ M Roper) + {^Roper/2y 

where Mn oper = 1440 MeV, T R oper = 325 MeV, and c is another normalization parameter. 



D. Compilation of the basic formulae 



This subsection contains a compilation of the main formulae used in the calculations and fits to the data. 
The elastic cross section is: 

da el . „,,.^^ 2(Q(t) - 1) 



at \ so 

The single diffraction (SD) dissociation cross section is 



A„F r \t) _ . (9) 



Double diffraction (DD) dissociation cross section: 



j3 / \ 2(a(t)-l) 

" -^72 - N DD F mel 2 (t,M 2 )F mel 2 (t,Ml) 



dtdMfdM* " m v ' 1 y v ' " ' 

with the norm iVo/j = 3-3—, with the inelastic vertex: 

F me 2 (t,M 2 ) = A res j^a^, p (Mf , t) + Cb g (TBg, 

where the Pomeron-proton total cross section is the sum N* resonances (Eq. (J7J) and the Roper resonance (Eq. 
with a relevant norm factor R (we remove the t dependent f re a(t) out of the sum): 

P 2 [fres(t)] 2 ■ M Roper ^ESspnJ ^ Xm a{M 2 ) 

v T (M.,t) = R ^ _ ^ + ^ v + [f reM n ^ 3(2n + . 5 _^ a(M|))2 V (Xma(M|))25 



and the background corresponding to non-resonance contributions: 



NB: /n Eq. [1$ M X) m p and m p are in [GeV\. 
The Pomeron trajectory is Q: 

a{t) = 1.075 + 0.34i, 
and the t— dependente elastic and inelastic form factors are: 

F p (t) = e b - 4 , f res (t) = e^"\ f b g(t) = e h »°\ 
The slope of the cone is defined as: 

where ^ stands for ^ff , ^p, or ^mc, defined by Eq. ©, (JIDJ and dHJ), respectively 

The local slope Bm at fixed M 2 is defined in the same way: 

d , 
= — In ■ 



tit ~ df dM 2 

NB: ^ in % {73^ is in units of [mb-GeV 2 ], and d f d ^,fi in Eq. U6]) in units of [mb-GeV 4 ]. 
The integrated cross sections are calculated as: 

do-so _ [ M * d 2 a SD ^ r2 

for the case of SD and: 



/ , ir o-dM 2 
dt J M 2 dtdMl 



^££ = / / f^D dM 2 xl dM 2 X2 



for the case of DD. 

We also calculate the fully integrated cross sections: 



-1 /-0.05s »2_ 

*sd= I dt / dM? dusD 



dtdMl 1 
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=J?JLj< dM ™^w (20) 



odd 

• / Ar)>3 - uiuivi- lU ivi- 2 

and 

d 2 a D D f 1 d 3 a D D 



dMldMl J dtdMldMl 



dt. (21) 
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V. RESULTS 



A. Fitting procedure 



The model contains 12 parameters, a large part of which is fixed either by their standard values (e.g. those of Regge 
trajectories, except for the Pomeron slope, whose slope exceeds the "standard" value to meet the SD data) or are set 
close to the previous fits [l|. 



• In our strategy we first adjust the model to the "standard candles" of elastic pp scattering at high energies 
(starting form 500 GeV). We considered only the data corresponding to the first cone, described by a linear 
exponential, implying also a linear Pomeron trajectory. 

The elastic data and Regge theory fix the parameters so,a(0), a', A e i, b e i. The relevant curves, the data and 
values of the fitted parameters are shown in Fig- EKa) , m Tables Hill and HT1 

The data at larger \t\, with the dip-bump structure and subsequent flattening of the cross section, both in elastic 
scattering and in SD may indicate the onset of new physics and the transition to hard scattering, implying a 
non-exponential residue and/or a non-linear Pomeron trajectory, see Refs. Q, that goes beyond the present 
study. 

NB: The parameter sq is strongly correlated with the slope parameters b e i, b res and bbg. 

• Single diffraction dissociation (SD) is an important pillar in our fitting procedure. The following parameters 
were fitted to SD data: A res , Cb g , R, b res , bb g , ?, r\. As input data we used: a)double differential cross sections 
^fP- V e rsus u 2 x at \t\ = 0.05 GeV 2 (see Fig. [TQa)) and b) at at \t\ = 0.5 GeV 2 (see Fig. HITb)): c)single 
differential cross sections dcr f t D vs. t (see Fig. [TOTb)); d) fully integrated cross sections versus energy y/s (see 
Fig-UDta)). 

In particular due to the difference between the kinematical regions chosen by the different experimental groups, 
there are some discrepancies int the data. 

At low t (below 0.5 GeV 2 ), the t— dependence of SD cross section are well described by an exponential fit, 
see Figs. ITOT b) andQTJb)), but beyond this region the cross sections start flattening due to transition effects 
towards hard physics. 

Low-energy data yfs < 100 GeV (see Figs. [TOTaL [To]) require the inclusion of non- leading Reggeons, so they are 
outside our single Pomeron exchange in the t channel. 

• Double diffraction dissociation (DD) cross sections follow, up to some fine-tuning of the parameters, 
from our fits to SD and factorization relation. 

Integration in Mf and Mf comprises the range A?y > 3, where A?/ = In I M s /?, a ) • 



The \ 2 values are quoted at relevant figures. The 
prediction are summarize in Table IVII 





Data [GeV -2 ] 


S ei (7TeV) 


19.9 ±0.3 [25] 


B e i(1.8 TeV) 


17.0 ±0.5 [26] 
17.9 ±2.5 [27] 
16.99 [28] 


B ei (546 GeV) 


15.35 [28] 
15.0 [29] 



of the fitted parameters are presented in Table IIVI and our 





Data [mb] 


Calcuation [mb] 


<M7TeV) 


25.4 ± 1.1 [25J 


24.5 


0^(1.8 TeV) 


16.6 ± 1.6 [26] 


17.99 


a ei (546 GeV) 


13.6 [30J 


13.8 



TABLE II: Forward slope of elastic pp scat- 
tering, see Fig. [9jb) 



TABLE III: pp elastic cross section, Eq. ([9]), calculated 
with the parameters quoted in Tab. IIVI 
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A e i [mb] 


33.579 


b el [GeV~ 2 ] 


1.937 


A res [mb-GeV 4 ] 


2.21 


C bg [mb] 


2.07 


R 


0.45 


6 res [GeV" 2 ] 


-0.507 


b bg [GoV- 2 ] 


-1.013 



so 


1 


c 


0.8 


n 


1 


a(t) = a(0) 


+ a't 


a(0) 


1.075 


Q'[GeV" 2 ] 


0.34 



TABLE IV: Fitted parameters, see Eqs. @, JTOj, (TXJ. 



°"Sd|5<0.05 


16.7 


a SD\ M<i GeV 


4.9 


a so\ M>i GeV 


11.8 


0"Dd|aj7>3 


8.6 


ctdd Im2a/|<i6 GeV 2 


0.4 


0-1315 Im?m?>i6 GeV 2 


8.2 



TeV] 


a SD \m 2 <200GsV 


O"S£>|{<0.08 


CDD Arj>3 


7 


13.7 


16.7 


8.6 


10 


14.6 


18.2 


9.4 


14 


15.5 


19.7 


10.2 



TABLE V: Cross sections [mb] calculated TABLE VI: Predictions (in [mb]) for SD and DD at 7, 10 and 
at 7 TeV. 14 TeV. 
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FIG. 9: (a) Elastic cross section (black line) caclulated from Eq. © and the TOTEM data points [25| at -\/s = 7 TeV; 
X 2 /ndf= 1.6. 

(b) Slope of the elastic cross section B = ^ (in (^f 1 )) • 
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C. Single diffraction dissociation 
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FIG. 10: (a) Single diffraction dissociation cross section vs. energy y/s calculated from Eq. (fl0|) . The points with red full 
circles correspond to (unpublished) the ALICE data (M < 200 GeV), the point with pink open circles are those extrapolated 
to ^ < 0.05 31]. For a\ (<0 . 05 : x7 n = 1-6, n= 6 (only y/1 > 100 GeV) and for <r| M<200 GeV : X 2 /n< 0.01, n= 3. 
(b) Single differential cross section a £ t D for SD calculated from Eq. (flOj) , integrated over the region Ml < 0.05s. 
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FIG. 11: Double differential cross sections for SD at t = —0.05 (a) and t - 
for 1800 GeV 2 and xV n = 4 -6, n= 8 for 546 GeV 2 (only Fig. (a)); 



-0.5 (b) calculated from Eq. ((TOj); xV n = 1-2, n= 8 
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FIG. 12: (a) Double differential SD cross sections as functions of M 2 for different t values, see Eq. (|10[) . 
(b) Double differential SD cross sections as functions of t for different M 2 values, see Eq. (|10[1 . 




M 2 (GeV 2 ) t (GeV 2 ) 

(a) (b) 

FIG. 13: (a) Single differential SD cross section as a function of M 2 integrated over t € [0.0 : 1.0], see Eq. (|10[1 . 

(b) Single differential SD cross section as a function of t, integrated over the M 2 from the threshold up to some fixed values, 

see Eq. JTOj) . 
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FIG. 14: (a) Slopes of differential SD cross section at different M 2 values, B = 4 In f °?3 



(b) Slopes of differential SD cross section integrated in M x over the region of the first resonance, then over the whole diffraction 
region. B = f t In where = ^v dM l 
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Double differential SD cross section j.Jf," is calculated from Eq. (jTUJ) at yfs = 7 TeV. 
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FIG. 15: (a) Normalized cross sections to showing the the prominent role of the resonances over the background (cf. with the 
finite mass sum rules, Fig. 26 in Ref. @|). 



(b) Double differential cross section 



for SD at 546 and 1800 GeV at |t| = 0.05. The blue points correspond to diffracion 



reggion 0.01 < £ < 0.05. The flattening of the cross section above £ > 0.05 may indicate the contribution of the Pomeron 
exchange in the Pp total cross section. However, this region may not be diffractive any more. 



D. Double diffraction dissociation 
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FIG. 17: Double differential DD cross section as a function of M\ and M| integrated over t; see Eq. (|11[) . 
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FIG. 18: (a) Single differential cross section da ® t D as a function of t integrated over the region of resonances {MlMl < 16 GeV 2 ) 
and over the whole region of diffraction (A77 > 3), see Eq. |QTJ|. (b) The slope B = ^ In (^§ s -) is calculated from Eq. (fTTj) . 




FIG. 19: Triple differential DD cross sections as a functions of M? and M| for different t values: i G {-0.1; -0.2; -0.3} 

calculated from Eq. (|lip . 
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VI. APPENDIX: THE KINEMATIC AL FACTOR K(x B ,MZ) AND TURN-DOWN OF THE CROSS 

SECTIONS AS t -> 



By guge invariance the DIS SF, and consequently our Pp transicion amplitude should vanish as Q 2 when Q 2 — > 0. 
This property is built in in the formula below, but in the main text, for simplicity we have ignored it. In this Appendix 
we show its consequences for the measurables (cross sections and slopes). 



We remind that the cross sections Eqs. (ITU|) . (fTTj) contain inelastic verteces with the kinematical factor: 

F mel 2 (t,M 2 ) = A res K f (xB,M 2 )a^ p (M 2 ,t) + C bg a Bg , 

xb(I-xb) 2 



(22) 



where Kf(xs,M 2 ) 



2-2 \ 3/2 • 



Here Kj{xb,M 2 ) relates the nucleon structure function to the Pomeron-proton total crass section (for details see 
Ref. [l[), xb — M 2~^2_ t is the Bjorken variable, that for high M 2 and moderate t is x s» —t/M 2 , and we remind 

that —t replaces the virtuality Q in proton-Pomeron DIS, and m p is the proton mass. For large s and small/moderate 
t, x fa —t/s. The Bjorken xb should not be confused with the Feynman variable xf, defined by 1 — xf ~ (M 2 — m 2 )/x. 

Below we repeat some results if the main text, but with the above kinematical factor included. 

With the kinematical factor included, some of the parameters should be refitted, namely: 

A res = 15.5 mb- GeV 2 , b res = -0.41 GeV~ 2 , b bg = -0.91 GeV" 2 . 

Other parameters are kept unchanged (see Tab. IIV|) . 
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FIG. 20: (a) SD cross section vs. energy y/s calculated from Eq. (|10p with the kinematical factor Eq. (|22[l . For cr|g<o.o5: 
X 2 /n= 1.6, n= 6 (only y/i > 100 GeV); for a| M<200 Ge y: X 2 /n= 0.02, n= 3. 

(b) DD cross section vs. energy y/s calculated from E q. (1111) with the kinematical factor Eq. (|22[) ; X 2 / n ~ 0.08, n — 7, for 
\fs > 100 GeV only. The experimental points are from [3ll ]. 
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FIG. 21: (a) Double differential SD cross sections as functions of M 2 for different t values, 
(b) Double differential SD cross sections as functions of t for different M 2 values. 
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FIG. 22: (a) Slopes of differential SD cross section at different M 2 values, B = 4 In f"?3 



Mj={Af2,M|,Af|} 



(b) Slopes of differential SD cross section integrated in M x over the region of the first resonance, and over the whole diffraction 
region. B = 4 In where %P = f 

° at dt 5 dt J 1 



IM 2 dtdM 2 x 

Double differential SD cross section d ^3 i s calculated from Eq. (|10|) with the kinematical factor Eq. (1221) for ^/s = 7 TeV. 
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FIG. 23: (a) Single differential cross section drT " t D as a function of t integrated over the region of resonances 
(M? Mi < 16 GeV 2 ) and over the whole region of diffraction (At] > 3). See Eq. (fTTjl with the kinematical factor Eq. ([22]) ; 
(b) The slope B = j- t In (^§ B -) calculated from Eq. (JTTJ) with the kinematical factor Eq. ((22)) . 
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VII. CONCLUSIONS 



Below is a brief summary of our results: 

• At the LHC, in the diffraction cone region (t < 1 GeV 2 ) proton-proton scattering is dominated (over 95%) by 
Pomeron exchange (quantified in Ref. [3]). This enables full use of factorized Regge-pole models. Contributions 
from non-leading (secondary) trajectories can (and should be) included in the extension of the model to low 
energies, e.g. below those of the SPS. 

• Unlike to the most of the approaches which use the triple Regge limit for construction of inclusive diffraction, our 
approaches based on the assumed similarity between the Pomeron-proton and virtual photon-proton scattering. 
The proton structure function (SF) probed by the Pomeron is the central object of our studies. This SF, similar 
to the DIS SF, is exhibits direct-channel (i.e. missing mass, M) resonances transformed in resonances in single- 
double- and central diffraction dissociation. The high-M behaviour of the SF (or Pomeron-proton cross section) 
is Regge-behaved and contains two components: one decreasing roughly like M~ m , m s» 2 due to the exchange 
of a secondary Reggeon (not to be confused with the Pomeron exchange in the t channel!). The latter dominates 
the large-M part of the cross sections. Its possible manifestation may be seen in the data see Fig. [IST b). On 
the other hand, the large-M region is the border of diffraction, £ > 0.05. 

• An important and intriguing prediction of the present model is the possible turn-down of the cross sections 
towards t = 0, see Appendix. The forward direction cannot be reached kinematically in SD or DD, moreover even 
the non-zero but small \t\ events are difficult to be reached, especially that they are masked by electromagnetic 
interactions (although weaker than in elastic scattering, see Ref. Further studies, both theoretical and 
experimental, of this intriguing phenomena are of great importance. 

• The results of our calculations, that are mainly predictions for the LHC energies 7, 8 and 14 TeV, are collected 
in Table fVTl Sec. FVl . Predicted values for integrated, within various limits of t and/or M 2 , cross sections can be 
also found in that section. The quality of the fit is quantified by the relevant x 2 values. 

• The model has an important and interesting prediction, following from the gauge invariance of the structure 
functions {Pp production amplitudes), namely that the cross sections turn down at very small values of \t\, 
probably accessible in the nearly forward direction of future measurements. This result was anticipated in 
Ref. p). 

• Our approach in this paper is inclusive, ignoring e.g. the angular distribution of the produced particles from 
decaying resonances. All resonances, except Roper, lie on the N* trajectory. Any complete study of the final 
states should included also spin degrees of freedom, ignored in the present model. 

• For simplicity we used linear Regge trajectories and exponential residue functions, thus limiting the applicability 
of our model to low and intermediate values of \t\. Its extension to larger \t\ is straightforward and promising. It 
may reveal new phenomena, such as the the possible dip-bump structure is SD and DD as well as the transition 
to hard scattering at large momenta transfers, although it should be remembered that diffraction (coherence) 
is limited (independently) both by t and £. 
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